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Abstract 



The paper presents a new theory of unfolding of eigenvalue surfaces of real 
Qh' symmetric and Hermitian matrices due to an arbitrary complex perturbation 

I near a diabolic point. General asymptotic formulae describing deformations of 

a conical surface for different kinds of perturbing matrices are derived. As a 
physical application, singularities of the surfaces of refractive indices in crystal 
optics are studied. 



1 Introduction 



Since the papers by Von Neumann and Wigner (1929)] and Teller (1937)| it is known 



that the energy surfaces in quantum physics may cross forming two sheets of a double 

cone: a diabolo. The apex of the cone is called a diabolic point, see Berry and Wilkinson (1984)| . 



This kind of crossing is typical for systems described by real symmetric Hamiltoni- 
ans with at least two parameters and Hermitian Hamiltonians depending on three 
or more parameters. From mathematical point of view the energy surfaces are de- 
scribed by eigenvalues of real symmetric or Hermitian operators dependent on pa- 
rameters, and the diabolic point is a point of a double eigenvalue with two linearly 
independent eigenvectors. In modern problems of quantum physics, crystal optics. 
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physical chemistry, acoustics and mechanics it is important to know how the diabohc 
point bifurcates under arbitrary complex perturbations forming topological singular- 
ities of eigenvalue surfaces hke a double coffee filter with two exceptional points or 
a diabolic circle of exceptional points, see e.g. Mondragon and Hernandez (1993)| , 
Mondragon and Hernandez (1996)| , [Shuvalov and Scott (2000)] , |Keck et al. (2003)| , 



Berry and Dennis (2003)] , [Korsch and Mossman (2003)| , [Berry (2004)] . 



In our preceding companion paper Mailybaev et al. (2004)| , a general theory of 
coupling of eigenvalues for complex matrices of arbitrary dimension smoothly depend- 
ing on multiple real parameters was presented. Two kinds of important singularities 
were mathematically classified: the diabolic points (DPs) and the exceptional points 
(EPs). DP is a point where the eigenvalues coalesce, while corresponding eigenvec- 
tors remain different (linearly independent), and EP is a point where both eigen- 
values and eigenvectors coalesce forming a Jordan block. General formulae describ- 
ing coupling and decoupling of eigenvalues, crossing and avoided crossing of eigen- 
value surfaces were derived. Both the DP and EP cases are interesting in applica- 
tions and were observed in experiments, see Ramachandran and Ramaseshan (1961)| , 
Dembowsky et al. (200 1)| , [Dembowsky et al. (2003)] , [Stehmann et al. (2004)| . 



In the present paper following the theory developed in Mailybaev et al. (2004)| 
we study effects of complex perturbations in multiparameter families of real symmet- 
ric and Hermitian matrices. In case of real symmetric matrices we study unfolding 
of eigenvalue surfaces near a diabolic point under real and complex perturbations. 
Origination of singularities like a "double coffee filter" and a "diabolic circle" is an- 
alytically described. Unfolding of a diabolic point of a Hermitian matrix under an 
arbitrary complex perturbation is analytically treated. We emphasize that the un- 
folding of eigenvalue surfaces is described qualitatively as well as quantitatively by 
using only the information at the diabolic point, including eigenvalues, eigenvectors, 
and derivatives of the system matrix taken at the diabolic point. 

As a physical application, singularities of the surfaces of refractive indices in crystal 
optics are studied. Asymptotic formulae for the metamorphoses of these surfaces 
depending on properties of a crystal are established and discussed in detail. Singular 
axes for general crystals with weak absorption and chirality are found. A new explicit 
condition distinguishing the absorbtion-dominated and chirality-dominated crystals 
is established in terms of components of the inverse dielectric tensor. 
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2 Asymptotic expressions for eigenvalues near a 
diabolic point 

Let us consider the eigenvalue problem 

Au = Au (1) 

for an m X m Hermitian matrix A, where A is an eigenvalue and u is an eigenvec- 
tor. Such eigenvalue problems arise in non-dissipative physics with and without time 
reversal symmetry. Real symmetric and complex Hermitian matrices correspond to 
these two cases, respectively. We assume that the matrix A smoothly depends on a 
vector of n real parameters p = (pi, . . . ,Pn)- Let Aq be a double eigenvalue of the ma- 
trix Aq = A(po) for some vector po. Since Aq is a Hermitian matrix, the eigenvalue 
Ao is real and possesses two eigenvectors ui and U2. Thus, the point of eigenvalue 
coupling for Hermitian matrices is diabolic. We choose the eigenvectors satisfying the 
normalization conditions 

(Ui, Ui) = (U2, U2) = 1, (ui,U2)=0, (2) 

where the standard inner product of complex vectors is given by (u, v) = Y^^i UiVi. 

Under perturbation of parameters p = Po + Ap, the bifurcation of Aq into two 
simple eigenvalues A+ and A_ occurs. The asymptotic formula for A± under multipa- 
rameter perturbation is Mailybaev et al. (2004)| 



^ ^ , (fll+f22,Ap) /(fii -f22,Ap)2 , ,n . wp A \ 

A± = Ao + ±y + (fi2, Ap)(f2i, Ap). (3) 

Components of the vector fjj = {f^j, . . . , f^j) are 

where the derivative is taken at po, and inner products of vectors in (jH)) are given by 
(a, b) = In expression (jHl) the higher order terms o(||Ap||) and o(||Ap||2) 

are neglected before and under the square root. Since the matrix A is Hermitian, 
the vectors fn and f22 are real and the vectors fi2 = f2i are complex conjugate. In 
case of real symmetric matrices A = A-'", the vectors fn, f22, and fi2 = f2i are real. 
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The asymptotic expression for the eigenvectors corresponding to X± takes the form 
Mailybaev et al. (2004)] 



a± fi2,Ap A±-Ao- f22,Ap 

U± = a±Ui + /?±U2, — = r 77 T-^ = 77 — . 5 

f3± A± - Ao - (fii, Ap) (f2i, Ap) 

Expressions (0) provide zero order terms for the eigenvectors u± under perturbation 
of the parameter vector. 

Now, consider an arbitrary complex perturbation of the matrix family A(p) + 
AA(p). Such perturbations appear due to non-conservative effects breaking symmetry 
of the initial system. We assume that the size of perturbation AA(p) ~ e is small, 
where e = ||AA(po)|| is the Frobenius norm of the perturbation at the diabolic point. 
Behavior of the eigenvalues A± for small Ap and small e is described by the following 
asymptotic formula Mailybaev et al. (2004)| 



(fll + f22, Ap) £11 + £22 

A± — Ao H \ 



, / ((fll ~ ^227 Ap) + £11 — £22)^ , f,r A \ , A \ , N 

± Y ^ + ((fl2, Ap) +£i2)((f21, Ap) +£21)- 



(6) 



The quantities Eij are small complex numbers of order e given by the expression 

= (AA(po)ui,Uj) . (7) 

A small variation of the matrix family leads to the following correction of the asymp- 
totic expression for the eigenvectors 

_ e , ne "± _ (fi2, Ap) +£12 _ A± - Ao - (f22, Ap) - £22 

u^-«^u, + /5^U2, ^-A^_Ao-(fn,Ap)-.n" (f2i,Ap) + .2i ' 

(8) 

The ratios a;^//5^ = a^_/l3^_ at the point of coincident eigenvalues A_|_ = A_. Hence, 
the eigenvectors u_,_ = u_ coincide, and the point of eigenvalue coupling of the per- 
turbed system becomes exceptional (EP). For some specific perturbations AA(p), the 
coupling point may remain diabolic under the conditions 

(fl2,Ap)+£i2 = 0, (f2l,Ap)+£21 = 0, (fn-f22,Ap)+£n-£22 = 0, (9) 

when both ratios in (jSJ become undetermined. 

We observe that asymptotic description of unfolding of diabolic singularity due 
to perturbation of the matrix family requires only the value of AA(p) taken at the 
coupling point po- Dependence of the perturbation AA on the vector of parameters 
p near the point po is not so important, since it influences higher order terms. 
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3 Unfolding of a diabolic singularity for real sym- 
metric matrices 



Let us assume that A(p) is an n-parameter family of real symmetric matrices. Then 
its eigenvalues A are real. Let Aq be a double eigenvalue of the matrix Aq = A(po) with 
two real eigenvectors Ui and U2. Under perturbation of parameters p = Po + Ap, the 
eigenvalue Aq splits into two simple eigenvalues A+ and A_. The asymptotic formula 
for A-t under multiparameter perturbation is given by equations and (0)), where 
the vectors fn, f22, and fi2 = f2i are real. Then, equation ^ takes the form 

(,,_,„_ «. + fe.Ap) y_ (f„-fe.Ap)^ ^^_^^^^^.^^ (10) 

Equation (fTUj) describes a surface in the space (pi,P2, • • • ,Pn, A), which consists of 
two sheets A+(p) and A_(p). These sheets are connected at the points satisfying the 
equations 

A± = Ao + ^(fii + f22,Ap), (fn - f22, Ap) = 0, (fi2,Ap) = 0, (11) 

where the eigenvalues couple: A+ = A_. Equations (jllll define a plane of dimension 
n — 2. Thus, the double eigenvalue is a phenomenon of codimension 2 in an n- 
parameter family of real symmetric matrices |Von Neumann and Wigner (1929)] ■ 

For the two-parameter matrix A(pi,p2) equation (fTUI) defines a double cone with 
apex at the point (po, Aq) in the space (pi,P2,A), see Figured The point (po, Aq) 
is referred to as a "diabolic point" Berry and Wilkinson (1984) due to the conical 
shape of the children's toy "diabolo". 

Let us consider a perturbation A(p) + AA(p) of the real symmetric family A(p) 
in the vicinity of the diabolic point po, where AA(p) is a complex matrix with the 
small norm e = ||AA(po)||. Splitting of the double eigenvalue Aq due to a change of 
the vector of parameters Ap and a small complex perturbation AA is described by 
equation (jHl), which acquires the form 

x± = \'o + fi±V^, c={x + 0' + {y + v)'-C- (12) 

In equation (fT^ the quantities Ag, x, and y are real: 

A[, = Ao + ^(fii + f22,Ap), a; = ^(fn-f22,Ap), y=(fi2,Ap), (13) 




Figure 1: A diabolic point in a family of real symmetric matrices. 

while the small coefficients fi, ^, rj, and ( are complex: 

1111 
= 2(^11 +^22), ^=2(^11-^22), ^? = 2(^12 + ^21), (=2(^12-^21). (14) 

Separating real and imaginary parts in equation (fT^ . we find 

Re^(A-A[,-;u)-Im^(A-A'o-/i) = Rec, 2Re(A-A'o-/i)Im(A-A'o-;u) = Imc, (15) 
where 

Rec= {lm\-lm^^-lm^ri-Re\)+{x+Re^y+{y+Reriy, (16) 

Imc = 2((x+ReOlm^ + {y+Rer])lm'r] - ReCImC). (17) 

From equations (fTH|) we get the expressions determining the real and imaginary 
parts of the perturbed eigenvalues 

ReA± = Ao + Re/i ± ^(^Rec + v^Re^c + Irn^c ^ /2, (18) 

ImA± = Im// ± ^ (^-Rec + \/Re^c + Wc ^ /2. (19) 

Strictly speaking, for the same eigenvalue one should take equal or opposite signs 
before the square roots in (jl8|) . (jl9p for positive or negative Imc, respectively. 

Equations (fTHj) and (fT^ define surfaces in the spaces (pi,p2, • • • ,Pn, ReA) and 
(pi,P2, • • • ,Pn, ImA). Two sheets of the surface (fTK|) are connected (ReA+ = ReA_) at 
the points satisfying the conditions 

Rec < 0, Imc = 0, ReA± = Aq + Re/x, (20) 
6 



while the sheets ImA+(p) and ImA_(p) are glued at the set of points satisfying 

Rec > 0, Imc = 0, ImA± = Im/i. (21) 

Note that in the neighborhood of the intersections ()2()|1 and 1)211] the eigenvalue 
sheets given by the formulae (fTHj) and (fT^ can be described by the following approx- 
imate expressions 



Imc / — 1 Imc / 1 

ReA±=An+Re/i ± — \/ — , Rec < 0; ImA±=Im/i ± — \/— , Rec > 0. (22) 
2 V Rec 2 V Rec 

The eigenvalue remains double under the perturbation of parameters when c = 0, 
which yields two equations Rec = and Imc = 0. Two cases are distinguished 
according to the sign of the quantity 

D = lm^^+Wr]-WC. (23) 

li D > 0, then the equations Rec = and Imc = with expressions (fTH|) . (fT7|) yield 
two solutions {xa,ya) and {xb,yb), where 



Im^ReCImC ± Im?7./(Im^^+Im^r7+Re^C)(Im^^+Im^r/-Im^C) 

Xa,b = T 2.,T 2 Re^, (24) 



Imr^ReCImC =F Im^V(Im^^+Wr/+Re^C)(Im^^+Wr/-WC) 

ya,b = T , T 2 ^^V- (25) 

These two solutions determine the points in parameter space, where double eigenvalues 
appear. When D = 0, the two solutions coincide. For D < 0, the equations Rec = 
and Imc = have no real solutions. In the latter case, the eigenvalues A+ and A_ 
separate for all Ap. 

Note that the quantities Im^ and Imrj are expressed by means of the anti-Hermitian 



part AAat = (AA - A A ) /2 of the matrix A A as 

(AAjv(po)ui, Ui)-(AAAr(po)u2, U2) 



Im^ 
Im?7 



2z 

[AA7v(po)ui, U2) + (AAAr(po)u2, Ui) 



(26) 



2i 



while lm( depends on the Hermitian part A A// = (AA + A A )/2 as 

^ (AAh(po)ui, U2) - (AAh(po)u2, ui) ^2^^ 
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Figure 2: Unfolding of a diabolic point due to complex perturbation. 

If -D > 0, one can say that the influence of the anti-Hermitian part of the perturbation 
A A is stronger than that of the Hermitian part. If the Hermitian part prevails in 
the perturbation AA, we have D < 0. In particular, D = —lm^(^ < for a purely 
Hermitian perturbation AA. 

Let us assume that the vector p consists of only two components pi and p2, and 
consider the surfaces (|18|) and (jl9|) for different kinds of the perturbation AA(p). 
Consider first the case D < 0. Then, the eigensheets ReA+(p) and ReA„(p) are 
separate, see Figure EK- Indeed, for D < — Re^^ the inequality Rec > holds for all 
variations of parameters, see equation (fT^ . In the case when — Re^C < D < the 
equation Rec = with expressions (fT!?|l define an ellipse in the plane of parameters 
(pi,P2)- Inside the ellipse we have Rec < and outside Rec > 0. Equation Imc = 
defines a line in parameter plane. The line and the ellipse have no common points 
for D < since there are no real solutions of the equation c = 0. Hence, for D < 
conditions (j^UI) are not fulfilled and the real parts of the eigenvalues avoid crossing. 
As the size of the complex perturbation decreases — 0), the two sheets come closer 
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and touch each other at the point (po, Aq) for e = forming the diabohc singularity. 
The sheets ImA+(p) and ImA_(p) of the eigensurface intersect along the line 

Imc/2 = (x+Re^)Im^ + {y+Reri)lmri — Re(lm( = 0, ImA± = Im/x, (28) 

given by conditions (PTj) . Note that, by using (j^ . one can show that the angle of 
intersection of the imaginary eigensheets is small of order e and tends to zero as e — 0. 

In the case D > the line Imc = and the ellipse Rec = have common points 
Pa and pb where the eigenvalues couple. Coordinates of these points found from the 
equations (|T!^ are 



where Xa,b and i/a^h are defined by expressions and ((23) • Here we have assumed 
that the vectors fn — ¥22 and fi2 are linearly independent. Note that the points pa 
and pb coincide in the degenerate case D = 0. 

According to conditions the real eigensheets ReA+(p) and ReA_(p) are glued 
in the interval [pa, p^] of the line 



Imc/2 = (x+ReOlm^ + {y+Reri)lmri - ReCImC = 0, ReA± = A[) + Re/z. (30) 



The unfolding of a diabolic point into the double coffee filter is shown in Figure |2t>. 

From conditions ()21|) it follows that the imaginary eigensheets ImA+(p) and ImA_(p) 
are connected along the straight line where the interval [pa, Pb] is excluded, see 
Figure |2b- According to the formulae ()22|) the angle of intersection of the imag- 
inary eigensheets tends to vr as the points pa and pb are approached, since Rec 
goes to zero. At far distances from the interval [p^, p^] this angle becomes small 
of order e. With the decrease of the size of complex perturbation e the interval 
shrinks and the angle of intersection goes to zero. At e = the imaginary parts 
of the eigenvalues coincide: ImA+ = ImA_ = 0. Note that in crystal optics and 
acoustics the interval [pa,Pf)] is referred to as a "branch cut", and the points Pa, Pb 
are called "singular axes", see IShuvalov and Scott (2000)| [Berry and Dennis (2003)| 
Ramachandran and Ramaseshan (1961)| . According to equation (jH]) the double eigen- 
values at Pa and pb possess only one eigenvector and, hence, they are exceptional 
points (EPs). 




/l2(/ll 122) fl2i.fll fh) /l2(/ll 122) /l2(/ll /2^2) 




(29) 



The surface of real eigenvalues (fTHj) is called a " double coffee filter" Keck et al. (2003)] . 
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Figure 3: A real non-symmetric perturbation of a diabolic point. 

Now let us consider the case when the perturbation AA(p) is real. In this case /x, 
^, T], (, and hence, 

c={x + 0' + {y + vf-C' (31) 

are real quantities. According to (fT^ the eigenvalues X± are complex-conjugate if 
c < and real if c > 0. The eigenvalues couple for c = forming a set consisting of 
exceptional points with double real eigenvalues. 

Consider a system depending on a vector of two parameters p = (pi,P2)- Then 
equation c = with expressions (fT^ and define an ellipse in parameter plane; 
c < inside the ellipse and c > outside. Real parts of the eigenvalues are given by 
the equations 

c > : (ReA-A[, - fif-{x + + vf= - (32) 

c < : ReA = Ao + /i. (33) 

Equation (jH^ defines a hyperboloid in the space (pi,P2, ReA). Real parts of the eigen- 
values A± coincide at the disk determined by equation (j^ . see Figure El Imaginary 
parts of the eigenvalues are 

c > : ImA = 0, (34) 
c<0: lm^X + {x + ^f+{y + rif = (35) 

The imaginary parts are both zero at the points of the plane ()34|1 surrounding the 
ellipsoid (jHS)) ("a bubble") in the space {pi,p2,lin.X), see Figure 01 The eigenvalues 
couple at the points of the elliptic ring 

A± = A'o + /i, {x + o'Hy + vy = C, (36) 
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consisting of exceptional points, see Figure El By that reason we call it an "excep- 
tional ring", which is a better name compared with a "diabolic circle" suggested by 



Finally, it is instructive to consider deformations of the surfaces fj32|l . (j33|l and 
()H4|1 . as the real perturbation becomes complex. If the imaginary part of the 
perturbation ImAA is such that D < 0, then the parts of the hyperboloid (jH^ 
connected by the disk are separated into the two smooth surfaces described by 
the equation (fTH|) . On the other hand, the ellipsoid surrounded by the plane 
is foliated into two sheets crossing each other along the line Imc = 0, see Figure 2a. 
Recall that the line Imc = does not intersect the ellipse Rec = 0. When D > 0, 
the disk foliates into two sheets crossing along the interval [pa,Pb], where the 
points Pa and pf, are given by expression ()29p. As the size of the imaginary part of 
the perturbation ImAA increases, the angle of intersection of real eigensheets grows. 
By this way, the purely imaginary perturbation deforms the hyperboloid ()32j) into the 
double coffee filter (fTH)) . see Figure Eb- The ellipsoid (jHSl) surrounded by the plane 
flH4|l is transformed into two smooth sheets intersecting along the line Imc = 0, where 
the interval [p^, p;,] is excluded. The angle of intersection grows as the size of the 
perturbation ImAA increases. 

4 Unfolding of a diabolic singularity for Hermitian 
matrices 

Let us consider a multi-parameter Hermitian matrix A(p). Assume that po is a 
diabolic point, where the matrix Aq = A(po) has a double real eigenvalue Aq with 
two eigenvectors. The splitting of Aq into a pair of simple real eigenvalues A+ and 
A_ is described by expressions ©,(01), where the vectors fn and £22 are real and the 
vectors fi2 = f2i are complex conjugate. By using expression 0, we find 



Mondragon and Hernandez (1993)] , [Mondragon and Hernandez (1996)] . 




(37) 
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where Ag, x, y, and z are real quantities depending linearly on the perturbation of 
parameters Ap as follows 

A[, = Ao + , = ifll^|^, ,= (Ref,„Ap), .= (lmf,„Ap). 

(38) 

The eigenvalues coincide iix = y = z = 0. Thus, if the system depends on three 
parameters and the real vectors fn — ¥22, Refi2, and lmfi2 are linearly independent, 
the eigenvalues A+ and A_ split for any nonzero perturbation Ap. For more than 
three parameters, the equations x = y = z = with relations ()38|) provide a plane 
of diabolic points in parameter space. This plane has dimension n — 3, which agrees 
with the well-known fact that the diabolic point is a codimension 3 phenomenon for 
Hermitian systems Von Neumann and Wigner (1929)j [Arnold (1972)| . 

Now let us consider a general non-Hermitian perturbation of the system A(p) + 
AA(p), assuming that the size of perturbation at the diabolic point e = ||AA(po)|| is 
small. The two eigenvalues A+ and A_, which become complex due to non-Hermitian 
perturbation, are given by asymptotic expressions (jHI), ©• With the use of the new 
coordinates (jHHj) . we write the expression (0) as 

X± = X'o + fx±V^, (39) 

where 

c=(x + 0' + (2/ + r/)' + (^-<)^ (40) 

and /i, ^, T], ( are small complex quantities of order e given by expressions (|T4|l . 
The eigenvalues couple (A+ = A_) if c = 0. This yields two equations 

Rec = (x + Re^Y + (y + Rerif + {z + ImCY - (W^ + Im^r^ + Re^C) = 0, (41) 

Imc = 2(lm^(x + ReO + lmT]{y + Ret]) - ReC(-2 + ImC)) = 0. (42) 
Equation (jlH) defines a sphere in (x, y, z) space with the center at (— Re^, —Re 77, — Im^) 



and the radius ylm^^ + Im^rj + Re^C, which are small of order e. Equation yields 
a plane passing through the center of the sphere. The sphere and the plane intersect 
along a circle. Points of this circle determine values of parameters, for which the 
eigenvalues A± coincide. Since c = at the coupling point, expression (jSI) for the 
eigenvectors takes the form 
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Figure 4: Unfolding of a diabolic point into an exceptional ring in parameter space. 

Thus, all points of the circle are exceptional points, where the two eigenvectors u_ 
and u+ merge in addition to the coupling of the eigenvalues A+ and A_. By using the 
linear expressions ()38|) . the set of exceptional points is found in the original parameter 
space p. The exceptional circle in (x, z) space is transformed into an exceptional 
elliptic ring in three-parameter space p, see Figure |3J 

Let us consider the plane (jl^ . at which the quantity c is real. By formula (jH^ . 
the real parts of the eigenvalues \± coincide inside the exceptional ring, where c < 0, 
and the imaginary parts of A± coincide outside the exceptional ring, where c > 0, see 
the dark and light shaded areas in Figure EJ 

We see that, under a general complex perturbation, a diabolic point of a three- 
parameter Hermitian system bifurcates into an exceptional ring. This ring has elliptic 
shape and grows proportionally to the size of perturbation e. The real and imaginary 
parts of the eigenvalues A± coincide, respectively, inside and outside the exceptional 
ring in the plane of the ring. 

Finally, let us study the stratification of parameter space given by the condition 
|Re (A+ — A_)| = const. For problems of quantum mechanics, this difference describes 
the size of a gap between two adjacent energy levels. By using expression (jHIIjl . we find 
(A+ — A_)^ = 4c. Separating real and imaginary parts in this equation and extracting 
Im (A+ — A_), we get 

Re^(A+ - A_) - 4Re2(A+ - A_)Rec - AWc = 0, (44) 

where Rec and Imc are given by the first equalities in ()4ip and ()42j) . Given fixed 
value of |Re (A+ - A_)|, equation (jSj) with dH}, and (jSHl) define an ellipsoid in 
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Figure 5: Surfaces corresponding to coincident real or imaginary parts of eigenvalues. 

three-parameter space enclosing the exceptional ring, see Figure Eti- Similar analysis 
provides the equation 

Im^(A+ - A_) + 4Im2(A+ - A_)Re c - Alm^c = (45) 

for a surface given by the condition |Im (A_|_ — A_)| = const. In three-parameter space 
equation (j^5|) defines a hyperboloid surrounded by the exceptional ring, see Figure Eb- 



5 Unfolding of optical singularities of birefringent 
crystals 

Optical properties of a non-magnetic dichroic chiral anisotropic crystal are character- 
ized by the inverse dielectric tensor 77, which relates the vectors of electric field E and 
the displacement D as Landau et al. (1984)| 



E = r/D. (46) 

A monochromatic plane wave of frequency u that propagates in a direction specified 
by a real unit vector s = (si, S2, S3) has the form 

D(r,t) = D(s)exp^u;f ^s^r-t) , (47) 
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where n(s) is a refractive index, and r is the real vector of spatial coordinates. With 
the wave ()47j) and the constitutive relation pHjl Maxwell's equations after some ele- 
mentary manipulations yield (see e.g. Berry and Dennis (2003)| ) 

1 



r7D(s) — ss r7D(s) 



77,2 fs 



-D s . 



(4^ 



Multiplying equation ()48|) by the vector S"^ from the left, we find that for plane 
waves the vector D is always orthogonal to the direction s, i.e., S"^D(s) = 0. By using 
this condition, we write (j48j) in the form of the eigenvalue problem 



fl 



u 



Au, 



(49) 



88-^)77(1 — SS"^) 

where A = 77"^, u = D, and I is the identity matrix. Since I — ss"^ is a singular matrix, 
one of the eigenvalues is always zero. Let us denote the other two eigenvalues by A+ 
and A_. These eigenvalues determine refractive indices n, and the corresponding 
eigenvectors yield polarizations. 

The inverse dielectric tensor is described by a complex non-Hermitian matrix 77 = 
Vtransp+ildichroic+'nchirai- The Symmetric part of -q consisting of the real matrix -q^^.^nsp 
and imaginary matrix T7rfjchroic constitute the anisotropy tensor, which describes the 
birefringence of the crystal. For a transparent crystal, the anisotropy tensor is real 
and is represented only by the matrix T/^ransp! f*^^ ^ crystal with linear dichroism it is 
complex. Choosing coordinate axes along the principal axes of ?7iransp5 '^^ have 



transp 



( T], 
V2 

V 



\ 


V3 J 



(50) 



The matrix 



"n dichroic 



I rif^ rif^ rif^ \ 
Vf2 V22 vis 

V Vfs ^23 vis J 



(51) 



describes linear dichroism (absorption). The matrix Tj^f^^^^i gives the antisymmetric 
part of T] describing chirality (optical activity) of the crystal. It is determined by the 
optical activity vector g = ((71, (72, fl's) depending linearly on s as 



'Hchiral 





93 

-92 



-93 


91 



92 

-91 




7s 



711 7i2 7i3 

712 722 723 
V 7l3 723 733 / 



S2 
\ S3 J 



(52) 
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where 7 is a symmetric optical activity tensor; this tensor has an imaginary part for 
a material with circular dichroism, see Berry and Dennis (2003)| for more details. 

First, consider a transparent non-chiral crystal, when ildichroic — 7 = 0. 
Then the matrix 

A(p) = (I - ss^)77,,„„,p(I - ss^) (53) 

is real symmetric and depends on a vector of two parameters p = (si, S2) (see 
Berry and Dennis (2003)] for other ways of introducing two parameters). The third 



component of the direction vector s is found as S3 = iy 1 — sf — S2, where the cases 
of two different signs should be considered separately. Below we assume that three 
dielectric constants f]i > f]2 > Vs different. This corresponds to biaxial anisotropic 
crystals. 

The nonzero eigenvalues X± of the matrix A(p) are found explicitly in the form 
[Lewin (1994)] 

tranp A I : 

(54) 



A- 



trace A _^ iy^2 trace (A^) - (trace A)^. 



The eigenvalues X± are the same for opposite directions s and — s. By using (j5(J|) and 
in (jS3)), it is straightforward to show that two eigenvalues A+ and A_ couple at 



So 



(S'i,S'2,S'3), Ao = r72; Si = ±^J (r/i - r/a)/ {% - r^a), 5*2 



0, 5, 



(55) 

which determine four diabolic points (for two signs of Si and S3), also called optic 
Ramachandran and Ramaseshan (1961)| . The double eigenvalue Aq = 772 of the 



axes 



matrix Aq = A(po), Po = (>S'i,0) possesses two eigenvectors 









/ 


^3 


\ 


Ui = 


1 


, U2 = 

















-Si 


J 



(56) 



satisfying normalization conditions Using expressions (jS^ and (|^. we evaluate 
the vectors fjj with components (jlj) for optic axes as 



fn = (0, 0), f22 = {2{V3 - Vi)Si, 0), fi2 = f2i = (0, - r^O^i^a). 



(57) 



By using ()55|) and (j37|) in (fTUI) . we obtain the local asymptotic expression for the cone 
singularities in the space (si, S2, A) as 



(A - r^2 - iV3 - Vi)Siisi - Si)y = (r/3 - r/i)^5f ((si - SiY + Sis^) 



{51 
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Figure 6: Diabolic singularities near optic axes and their local approximations. 

Equation is valid for each of the four optic axes 

As an example, consider the case of rji = 0.5, ri2 = 0.4, rj^ = 0.1. Conical surfaces 
fl58|) are shown in Figure IHl together with the exact eigenvalue surfaces ()54|) . The two 
optic axes presented in FigureElare Sq = (±1/2, 0, -\/3/2) with the double eigenvalue 
Ao = 2/5; the eigenvalue surfaces for the opposite directions Sq = (±1/2, 0, — -\/3/2) 
are exactly the same. 

Now let us assume that the crystal possesses absorption and chirality. Then the 
matrix family ()53|) takes a complex perturbation A(p) + AA(p), where 

AA(p) = (I - ss^){rjaichroic + Vchirai)(i - ss'^)- (59) 

Assume that the absorption and chirality are weak, i.e., e = HVdichroicW ± ll^c/iira«ll is 
small. Then we can use asymptotic formulae of Sections 2 and 3 to describe unfolding 
of diabolic singularities of the eigenvalue surfaces. For this purpose, we need to know 
only the value of the perturbation AA at the optic axes of the transparent non-chiral 
crystal Sq. 

Substituting matrix (jHH) evaluated at optic axes into expression (|7j), we obtain 

ei2 = + liiSi + inS;)S^ + i{r,i^ - Tis^i - 733^3)^3, (60) 
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By using formulae (HH), we get 

/i = z(r^2^2 + vfiSl - 2vf,S,Ss + vi3S!)/2, 
e = ^{vi2 - vfiSl + 2vf,SiS, - vi3S!)/2, 

(61) 

1^ 

C = -^ (7ll5? + 2713^153 + 733^1)- 

We see that fi, ^, and rj are purely imaginary numbers depending only on dichroic 
properties of the crystal (absorption). The quantity ( depends only on chiral proper- 
ties of the crystal; ( is purely imaginary if the optical activity tensor 7 is real. 

Singularities for crystals with weak dichroism and chirality were studied recently 
Berry and Dennis (2003)| . It was shown that the double coffee filter singularity 



m 



arises in absorption-dominated crystals, and the sheets of real parts of eigenvalues 
are separated in chirality-dominated crystals. According to the results of Section 3, 
these two cases are explicitly determined by the conditions D > and D < 0, re- 
spectively, where D = Im^^ -|- Im^r] — lia^C- These conditions are new and important 
because they provide quantitative definitions of absorption-dominated and chirality- 
dominated regimes for unfolding of the diabolic singularity in terms of components 
of the inverse dielectric tensor. Indeed, according to ()61|). ^ and 77 depend linearly on 
all the components of the tensor "naichroicy while ( depends linearly on the components 
lijy hj = 1, 3 of the optical activity tensor 7. 

Note that according to the sign of the quantity D taken at different optic axes we 
can classify crystals by their optic properties. For example, the important case is a 
chirality-dominated crystal with D < for all four optic axes. Then real parts of the 
eigenvalues separate for all directions s. 

There are four optic axes (jKKjl . which determine two pairs of opposite space direc- 
tion ±So. It is easy to see that the unfolding conditions coincide for the optic axes 
given by opposite directions, while these conditions are different for different pairs 
of optic axes. In the absorption-dominated case, when diabolic singularities unfold 
into coffee-filters near two opposite optic axes ±So = ±(»S'i, 0, 5*3), the four exceptional 
points of eigenvalue coupling ±Sa and ±85 (also called singular axes) appear. By using 
(|37j) in (fT^ . we obtain the asymptotic formulae 



sf = S,+ ^"'^ sf = - sf = ^^l-isff-isff, (62) 

(^1 - ri3)Si ivs - VijSiSs V 
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for the components of the vectors s^,;,, where Xa,b and ya,b are found by using ex- 
pressions fl24j) ■ ()25|1 . and (pTT| . In particular, for non-chiral crystals, we have C = 0. 
Then expressions (j^ . yield Xa,b = ilmr/ = r7f2'S'3 — '723*S'i and yafi = T^m^ = 



Equation Imc = determines a line of singularities in the parameter space p = 
(si, S2). By using (jFTf) . in (fT^ . (fTTj) . we find this line in the form 

(si - Si)Si{r]i - ?73)Im^ - S2SiS-i{r]i - 773) \mr] - ReCImC = 0. (63) 

In the absorption-dominated case, line (|FK?jl contains two exceptional points pa,6 = 
{si^, 82^) corresponding to the singular axes A segment between the points 
and pf, corresponds to the coincidence of real parts of the eigenvalues Re A+ = Re A_, 
while imaginary parts of the eigenvalues ImA+ = ImA_ merge at points of line ()63p 
outside this segment, see Figure |2|3. In the chirality-dominated case, when singular 
axes do not appear, imaginary parts of the eigenvalues ImA+ = ImA_ coincide at 
points of the whole line ()63|) . see Figure Et- If the optical activity tensor 7 is real or 
purely imaginary, then the line of singularities ()fi3j) passes through the diabolic point 
Po, and position of this line does not depend on 7. 

As a numerical example, let us consider a crystal possessing weak absorption and 
chirality described by the tensors (jSH), with 



(64) 



A corresponding transparent non-chiral crystal is characterized by r/i = 0.5, ri2 = 0.4, 
rjs = 0.1, and its eigenvalue surfaces with two optic axes are presented in Figure IHl 
By using (jMj) in (jHH), we find that the condition D = igQQoo (4v^ — 5) > is satisfied 
for the left optic axis Sq = {—1/2, 0, v^/2). Hence, the diabolic singularity bifurcates 
into a double coffee filter with two exceptional points whose coordinates according to 
expressions fjH^ are 
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(65) 



Local approximations of the eigenvalue surfaces are given by expressions (fTK|) . (fTIUl . 
where 



Re c 



35 - 28^3 
160000 



+ —(si + 1/2)2 + 

25^ ' ' 100 



'2' 



Imc 



6 + V3 
2000 



S2- 



(66) 
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Figure 7: Unfolding of singularities near optic axes. 

Figure UK shows these local approximations compared with the exact eigenvalue sur- 
faces given by (jMjl . For the right optic axis Sq = (1/2, 0, V3/2), the condit ion 
D = — j^gQQQQ (4v^ + 5) < is satisfied. Hence, the eigenvalue sheets (for real parts) 
separate under the bifurcation of the right diabolic singularity. Approximate and 
exact eigenvalue surfaces are shown in Figure Ub- The approximations are given by 
expressions f|THj) . (fT^ . where 

35 + 28^3 1, ,,^,2 3 2 T 6-^3 
Rec= (si- 1/2)2 H ^2 jj^g^ _g (57^ 

160000 25^ ' ^ 100 2 2000 ^ ^ 

We observe that the unfolding types are different for different optic axes. As it is seen 
from Figure [71 the asymptotic formulae provide an accurate description for unfolding 
of eigenvalue surfaces near diabolic points. 



6 Conclusion 

Non-Hermitian Hamiltonians and matrices usually appear in physics when dissipative 
and other non-conservative effects are taken into account. The known examples are 
complex refractive indices in optics and complex potentials describing the scattering 
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of electrons or X-rays. Traditionally, non-Hermitian matrices appear in physics as 
a perturbation of Hermitian matrices. As it is stated in Berry (2004)| , Hermitian 
physics differs radically from non-Hermitian physics in case of coalescence (coupling) 
of eigenvalues. In the present paper we have studied this important case carefully. 
We gave analytical description for unfolding of eigenvalue surfaces due to an arbitrary 
complex perturbation with the singularities known in the literature as a " double coffee- 
filter" and a " diabolic circle" . We emphasize that the developed theory requires only 
eigenvectors and derivatives of the matrices taken at the singular point, while the 
size of the matrix and its dependence on parameters are arbitrary. This makes the 
presented theory powerful and practical for a wide class of physical problems. The 
given physical example from crystal optics demonstrates applicability and accuracy 
of the theory. 
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